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摘 要 : 喀 喇 昆仑 山 分 布 有 众多 的 跃 动 型 冰川 , 跃 动 机 理 存在 异 质 性 。 使 用 Landsat Sentinel-1A、 


TSX/VTDX 等 多 源 遥 感 数据 , 获 取 了 中 喀 喇 昆仑 布 拉 尔 杜 冰川 在 跃 动 期 间 冰 川 表 
台 , 流 速 逐渐 增加 ;到 2013 年 之 后 进入 快速 运动 期 ,最 
速度 急剧 下 降 , 随 后 保持 较 低 的 流速 至 2016 年 1 月 ， 


变化 。 结 果 表 明 : 布 拉 尔 杜 冰 川 从 2006 年 
高 流速 可 达 4.9 md';2015 年 8 月 底 , 冰 川 表 鱼 


高 程 和 表面 流速 


而 后 流速 再 次 增加 ,到 同年 2 月 初 便 进入 平静 期 。2000 一 2014 年 冰川 主干 中 上 游 有 明显 隆起 ,而 冰 
川 接收 区 明显 减 薄 , 最 大 减 薄 达 89 m;2014 一 2018 年 冰川 主 于 中 上 游 以 及 各 支流 均 有 不 同 程度 的 


减 薄 ,冰川 主干 的 接收 区 高 程 显著 增加 ,最 大 增 厚 120 mo RAK II K E 


高 程 变 化 以 及 流速 变化 的 


特征 ,认为 布 拉 尔 杜 冰川 的 支流 引发 了 此 次 跃 动 , 且 本 次 跃 动 受 水 文 机 制 的 影响 较 大 ;结合 现 有 的 
数据 和 文献 ,推断 布 拉 尔 杜 冰川 的 跃 动 间隔 约 为 40a; 为 喀 喇 昆仑 冰川 跃 动 研究 提供 更 多 的 实例 ， 


也 可 为 此 区 域 冰 川 灾害 预警 研究 提供 参考 。 


X 键 词 : 冰川 跃 动 ; 表面 流速 ; RASE; TanDEM-X; 布 拉 尔 杜 冰川 


文章 编号 : 


冰川 牙 动 是 冰川 准 周期 性 的 运动 ,在 短暂 的 路 
动 阶段 (通常 以 高 流速 为 特征 ) 和 漫长 的 静止 阶段 
(以 末端 停滞 或 退缩 为 特征 ) 之 间 交 在" 于" EERS 
阶段 ,大 量 的 冰川 物质 从 上 游 积蓄 区 快速 向 下 游 接 
收 区 迁移 ,导致 冰川 表面 高 程 发 生 剧烈 变化 。 大 规 
模 的 冰川 物质 再 分 配 可 能 会 导致 冰川 末端 前 进 、 冰 
川 湖 演 决 洪水 (GLOF) 等 一 系列 灾害 ,从 而 威胁 流域 
内 居民 的 生命 和 财产 安全 。 例 如 ,2015 年 东 帕 米尔 
地 区 的 克拉 牙 依 拉克 冰川 发 生 跃 动 ,引发 的 雪崩 造 


山 、 帕 米尔 等 中 亚 部 分 山系 ”” ,俄罗斯 的 高 加 索 和 
新 西伯 利 亚 群 岛 等 地 区 也 曾 出 现 过 关于 冰川 跃 动 
的 报道 ”。 虽 然 跃 动 型 冰川 仅 约 占 全 球 冰 川 总 数 
的 1%, 但 是 其 重要 影响 却 不 可 和 忽视。 

目前 存在 两 种 普遍 接受 的 假说 解释 冰川 跃 动 : 
冰 下 水 文 机 制 模型 和 冰 下 热 转换 机 制 模型 冰 
下 水 文 机 制 模型 是 在 阿拉 斯 加 杂 色 (Variegated ) 冰 
JI] 1982—1983 年 跃 动 研 究 的 基础 上 发 展 起 来 的 ™， 
水 文 控制 的 跃 动 主要 由 于 冰 下 孔 际 水 压力 的 改变 


成 大 量 的 牧场 和 牧民 的 房屋 被 破坏 ”;2016 年 西藏 
阿 鲁 错 冰 川 跃 动 / 冰 衣 造成 牧民 和 牲畜 死亡 *。 
此 ,研究 冰川 跃 动 的 过 程 及 控制 机 理 对 于 预防 和 减 
轻 跃 动 带 来 的 灾害 具有 重要 理论 和 实践 意义 。 

跃 动 型 冰川 并 未 均匀 的 分 布 在 世界 冰川 区 , 它 
们 往往 聚集 在 少数 地 区 一 ,如 阿拉 斯 加 及 加 拿 大 育 
` 斯 瓦尔 巴 德 群岛 .格陵兰 岛 以 及 喀 喇 昆仑 .天 


Ht 
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而 导致 的 不 稳定 ,通常 伴随 着 冰 下 排水 通道 的 重 
组 ,这 类 跃 动 通常 在 排水 效率 低下 的 冬季 开始 SH 
束 于 排水 效率 高 的 夏季 。 热 转换 机 制 模型 是 在 斯 
瓦尔 巴 德 地 区 冰川 的 研究 中 发 展 起 来 的 ” , 热 转换 
机 制 的 跃 动 是 由 于 冰川 积蓄 区 的 物质 不 断 积累 致 
使 冰川 底部 温度 达到 融 点 ,冰川 底部 静水 压力 的 上 
升 致使 底部 快速 滑动 , 当 底部 融 水 全 部 流失 后 跃 动 
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结束 。 由 于 冰川 底 厂 特 性 很 难 观测 ,大 部 分 地 区 的 
跃 动机 制 仍 不 明确 。 

在 全 球 变 上 暧 的 气候 背景 下 ,全 球 冰川 整体 上 号 
现 不 同 程度 的 持续 退缩 , 而 喀 喇 昆仑 地 区 由 于 有 异常 
的 区 域 物质 平衡 “和 大 量 唉 动 型 冰川 “ 售 受 关注 。 
研究 表明 喀 喇 昆仑 的 跃 动机 理 存在 异 质 性 ,一 方面 
认为 咯 喇 昆仑 地 区 的 跃 动 是 受热 转换 机 制 控制 ， 
Quincey 等 "利用 卫星 遥感 对 喀 喇 昆仑 地 区 5 条 冰 
川 的 跃 动 过 程 进行 研究 ,量化 了 这 5 条 冰川 流速 的 
时 空 变化 ,并 认为 喀 喇 昆仑 地 区 近期 的 冰川 跃 动 是 
由 热 转 换 机 制 控 制 的 , 且 与 该 地 区 降水 增加 和 冰川 
积累 模式 导致 的 高 海拔 地 区 变 暖 相 一 致 ; 另 一 方面 
对 单个 冰川 路 动 的 事件 观测 表明 冰 下 水 文 条 件 的 
变化 是 主要 触发 机 制 , 如 Bhambri 等 ' 引 基于 多 源 遥 
感 数 据 对 喀 喇 昆仑 地 区 Shispare 冰川 跃 动 前 后 的 流 
速 .高 程 和 表面 特征 的 变化 进行 分 析 , 认 为 Shispare 
冰川 的 跃 劲 是 由 冰 下 水 文 机 制 控 制 。 

目前 喀 喇 昆仑 地 区 详细 的 冰川 跃 动 过 程 描述 
仍然 缺乏 ,至 今 仅 有 Shispare Hispar Kyagar, Khur- 
dopin 以 及 Baltoro 等 冰川 的 跃 动 过 程 有 相对 详细 的 
和 遥感 观测 “”。 随 着 遥感 技术 快速 的 发 展 使 得 对 跃 
动 冰 川 的 识别 和 研究 更 加 快捷 ,利用 获得 的 冰川 路 
动 过 程 及 特征 可 以 推断 冰川 的 雅 动 控制 机 制 。 本 
研究 收集 了 布 拉 尔 杜 冰 川 跃 动 期 间 Sentinel-1 数 


图 例 
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王 振 峰 等 : 中 喀 喇 昆仑 布 拉 尔 村 冰川 近期 跃 动 分 析 


据 ,基于 雷达 影像 偏 移 追 踪 法 提取 了 此 冰川 表面 运 
动 速度 ,此 外 还 基于 TerraSAR-X/TanDEM-X 数 据 生 
成 的 数字 高 程 模型 (DEM ) 对 冰川 表面 高 程 变化 进 
行 量化 ,描绘 其 跃 动 特征 ,进而 对 其 控制 机 理 进行 
探讨 ,以 期 对 喀 喇 昆仑 布 拉 尔 杜 冰川 跃 动 的 进一步 
认识 ,为 喀 喇 昆仑 冰川 跃 动 研究 提供 更 多 的 实例 。 


1 研究 区 概况 


布 拉 尔 杜 冰川 位 于 喀 喇 昆仑 中 部 ,其 融 水 汇 和 人 
什 玛 莎 河 , 上 游 东 侧 与 音 苏 盖 提 冰川 相 接 。 喀 喇 昆 
仑 地 区 主要 受 亚洲 季风 影响 ,夏季 降水 量 占 总 降水 
量 的 80%。 布 拉 尔 杜 冰川 编号 为 RGI60-14.04411， 
中 心 位 置 为 36%07'N,75°52'E( 图 1), 此 冰川 主体 呈 
现 自 南 向 北 的 流向 ,是 由 多 条 支 冰川 组 成 的 树枝 状 
山谷 冰川 ,面积 约 为 176 km””。 冰 川 末端 被 大 量 表 
磷 所 覆盖 ,表面 还 育 有 少量 冰 面 湖 。 


2 数据 与 方法 


2.1 数据 来 源 

2014 年 发 射 的 Sentinel-1A 卫星 搭载 C 波 段 传 
感 器 ,具有 全 球 采集 . 重 访 周期 短 和 良好 的 空间 分 
辩 率 等 特点 ,使 其 成 为 研究 冰川 运动 的 重要 数据 
源 。 本 研究 利用 Sentinel-1A 干 涉 宽 幅 (IW) 模 式 的 


注 :图 中 数字 为 支 冰川 编号 。 
图 1 布 拉 尔 杜 冰川 位 置 
Fig. 1 Location of the Braldu Glacier 
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影像 获取 布 拉 尔 杜 冰川 2014 一 2018 年 的 流速 变化 ， 
并 结合 ITS_LIVE 项 目 (https://nsidc.org/apps/itslive/) 
中 的 高 亚洲 地 区 冰川 的 年 均 流 速 数 据 进 行 分 析 。 

TerraSAR-X/TanDEM-X (TSX/TDX ) 系 统 是 由 德 
国 宇 航 局 (DLR ) 先 后 发 射 的 两 颗 X 波 段 SAR 卫星 组 
成 的 单轨 双 天 线 双 星 分 布 式 系统 。 两 颗 卫 星 具 有 
相似 的 结构 , 且 采 用 近 距 离 螺旋 结构 运行 ,重复 周 
期 为 11 4d。 该 系统 在 双 基 站 模式 下 ,可 以 同时 采集 
干涉 对 的 主 从 影像 ,具有 “0 时 间 基 线 ” 的 特点 。 
SRTM 是 在 2000 年 由 美国 国家 航空 航天 局 和 国家 影 
像 与 测绘 局 联合 进行 的 航天 飞机 雷达 地 形 测绘 任 
务 ,SRTM 的 观测 宽度 为 225 km, 其 生成 了 全 球 约 
80% 陆 地 表面 的 数字 高 程 模型 。 本 文 使 用 Ter- 
raSAR-X/TanDEM-X 影像 提取 新 的 DEM, = SRTM 
DEM 共同 用 于 计算 布 拉 尔 杜 冰川 跃 动 前 后 表面 高 
程 的 变化 ,同时 SRTM DEM 也 作为 参考 DEM 对 SAR 
数据 进行 基于 地 形 的 精 配 准 。 

Landsat 数 据 通过 美国 地 质 调 查 局 (USGS ,https: 
//earthexplorer.usgs.gov/) 进 行 处 理 和 分 发 。 本 文选 
取 了 夏秋 两 季 云 雾 较 少 的 影像 ,主要 用 于 监测 布 拉 
尔 杜 冰川 跃 动 期 间 冰 川 末端 和 内 侧 冰 研 的 变化 。 
本 研究 中 使 用 数据 见 表 1。 

2.2 冰川 表面 流速 与 高 程 提取 

2.2.1 冰川 表面 流速 提取 ”由 于 冰 流 速度 过 快 而 导 
致 的 冰 面 失 相 干 使 得 D-InSAR 难以 获得 高 质量 的 干 
涉 条 纹 图 ,而 偏 移 量 跟踪 法 中 的 强度 追踪 法 却 可 以 
克服 冰 面 失 相干 的 问题 ,被 认为 是 获取 冰川 流速 的 
可 靠 手 段 *””。 强 度 追 踪 的 原理 是 根据 图 像 的 纹理 
《条纹 或 后 向 散射 强度 ) 的 最 大 相关 系数 来 计算 目 
标的 偏 移 量 。 本 文采 用 瑞士 GAMMA 合成 孔径 雷达 
软件 中 的 Offset Tracking 模 块 对 Sentinel-1A 影像 进 
行 偏 移 量 估计 ,具体 流程 如 图 2。 在 偏 移 量 追踪 计 
算 时 需要 对 搜索 窗口 的 大 小 进行 设 定 ,窗口 尺寸 过 


小 会 导致 单个 窗口 不 稳定 ,可 能 无 法 在 其 他 图 像 中 
找到 对 应 的 窗口 。 经 试验 ,将 匹配 窗口 设置 为 
256x64 效 果 最 佳 。 

本 文 利用 SAR 影像 偏 移 量 追踪 技术 获取 的 冰 
川 流 速 ,其 误差 主要 来 自 图 像 配 准 图像 对 的 偏 移 
估计 及 系统 误差 。 为 了 对 冰川 表面 流速 的 精度 进 
行 评估 ,通常 假设 非 冰 川 区 为 静止 区 域 , 静 止 区 域 
包含 上 述 所 有 误差 源 。 因 此 通过 非 冰 川 区 残 差 位 
移 对 冰川 流速 结果 进行 评估 ,具体 公式 * 为 : 


e = VMED* + SE” (1) 
Wen 
N oa R 
Nu = ID (3) 


式 中 : eu 为 非 冰 川 区 偏 移 量 的 误差 ;MED 为 非 冰 川 
区 偏 移 量 的 平均 值 ;SE 为 非 冰 川 区 偏 移 量 的 标准 误 
差 ;SD 为 非 冰 川 区 偏 移 量 平均 值 的 标准 差 ; Ng 为 
去 自 相 关 后 的 像 元 个 数 ; Ni,, 为 非 冰 川 区 像 元 总 
数 ;及 为 像 元 分 辨 率 ;了 为 去 除 空间 自 相 关 距 离 。 根 
据 上 述 公 式 可 得 非 冰 川 区 偏 移 量 误差 为 0.016 med”, 
远 小 于 冰川 速度 的 观测 值 。 

2.2.2 冰川 表面 高 程 提取 及 变化 计算 在 瑞士 
GAMMA 软件 平台 下 对 TSXTDX SAR 数据 进行 差分 
干涉 处 理 , 在 差分 干涉 处 理 的 过 程 中 引入 1 rades” 
SRTM DEM 作为 参考 DEM, 从 而 获取 TSX/TDX 
DEM。 主 要 处 理 步骤 :首先 将 SRTM DEM 模拟 转换 
H SAR 成 像 数 据 , 其 次 TSX/TDX =j SRTM 模拟 的 
SAR 数据 进行 配 准 并 共 斩 相 乘 生 成 干涉 网 ” ,通过 
差分 处 理 , 去 除 SRTM DEM 模拟 的 地 形 相 位 和 基线 
误差 导致 的 线性 相位 ,从 而 得 到 冰川 表面 高 程 变化 
的 残 差 相位 ,将 得 到 的 残 差 相位 加 上 SRTM 模拟 的 
地 形 相 位 ,最 后 经 过 地 理 编码 , 即 可 生成 高 精度 的 
TSX/TDX DEM'”* ,具体 流程 如 图 3。 将 生成 的 TSX/ 


表 1 使 用 的 数据 源 


Tab.1 Datasets used in this study 


数据 用 途 分 辩 率 日 期 
Sentinel-1A 冰川 流速 监测 5 mx20 m 2014-10-31 一 2018-03-21 
TerraSAR-X/TanDEM-X 冰川 高 程 变 化 监测 1.4 mx2.2 m 2014-02-19 .2018-08-27 
SRTM-C 冰川 高 程 变化 监测 30 mx30 m 2000-02 
Landsat/OLI 冰川 表面 变化 分 析 30 mx30 m 2013-07-05 2014-07-24 ,2015-09-13 ,2016-10-01 „2017-09-20 
Landsat/ETM+ 冰川 表面 变化 分 析 30 mx30 m 2005-09-02 .2006-07-26 .2007-07-09 ,2012-07-10 
Landsat/TM 冰川 表面 变化 分 析 30 mx30 m 2008-09-25 „2009-08-27 „2010-10-17 ,2011-08-17 
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TDX DEM 4 SRTM DEM 进行 配 准 后 ,采用 大 地 测量 
法 计算 2000 一 2014 年 .2014 一 2018 年 布 拉 尔 杜 冰 川 
的 表面 高 程 变化 。 因 SRTM-X 未 完全 覆盖 研究 区 ， 
故 本 文选 取 SRTM-C 进行 研究 。 由 于 微波 对 冰雪 具 
有 一 定 的 穿 透 性 ,利用 InSAR 得 到 的 冰川 表面 高 程 
可 能 低 于 实际 高 程 。 考 虑 到 TSX/TDX 与 SRTM-X 的 
载 频 几 乎 相同 ,因此 在 进行 冰川 高 程 变 化 计算 时 ， 


SRTM DEM 


主 影像 
SLC 


、| 冰川 表 
面 流速 


注 :SRTM DEM 为 航天 飞机 雷达 地 形 测绘 获取 的 数字 高 程 模型 ; 
SLC 为 单 视 复 数 影像 。 下 同 。 
图 2 偏 移 量 追 踪 方 法 获取 冰川 表面 流速 流程 


Fig.2 Schematic diagram for acquirement of the surface 


辅 影像 | 
SLC 


B+ MR St 
TRES 


velocity of glaciers using offset tracking method 


SRTM DEM | 


TerraSAR-X 


注 :TSX/TDX 为 双 基 站 条 带 模式 获取 的 单 视 复 数 影 像 ;TSX/TDX 
DEM 为 TSX/TDX 经 差分 干涉 生成 的 数字 高 程 模型 。 
图 3 TSX/TDX 差 分 干涉 获取 DEM 流程 
Fig.3 Flow chart for DEM generation from TSX/TDX 


differential interferometry 


相 
位 
转 
[可 
各 
z 


250 F (a) 2000—20144F 
MED=0.84 
200 上 SE=0.11 
o =0.85 
150 + 
x x 
N EN 
100 上 
50 | | 
0 1 alll 川 h 1 1 
-20 -10 0 10 20 
高 程 变 化 /m 


注 :MED 为 非 冰 川 区 偏 移 量 的 平均 值 ;SE 为 非 
图 4 非 冰 川 区 表 
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可 利用 SRTM-X DEM 对 C 波 段 的 冰雪 穿 透 深度 进行 
估算 和 校正 ?2 ,经 计算 得 出 SRTM-C 在 喀 喇 昆仑 地 
区 的 穿 透 深度 约 为 2.4 m, 

冰川 表面 高 程 变化 的 不 确定 性 评估 方法 与 上 
文 冰川 表面 速度 的 不 确定 性 评估 类 似 , 假 定 非 冰 川 
区 高 程 无 变化 ,利用 非 冰 川 区 的 高 程 差 对 DEM 数据 
间 相 对 误差 进行 评估 。 最 终 不 同 DEM 数据 间 高 程 
误差 (c) 可 用 下 列 公 式 表示 52: 


o = MED?+SE? (4) 
SE= SD (5) 
JN 


式 中 :N 为 去 自 相 关 后 的 像 元 个 数 , 经 计算 得 到 
2000 一 2014 年 和 2014 一 2018 年 非 冰 川 区 表面 高 程 
谈 化 的 总 体 误差 分 别 为 0.85 m 和 2.41 m( 图 4)。 


3 结果 与 分 析 


3.1 冰川 表面 流速 变化 

利用 偏 移 量 追踪 的 方法 获取 了 布 拉 尔 杜 冰川 
2014 年 10 月 一 2018 年 1 月 的 时 序 表面 流速 数据 ,为 
了 进一步 地 观察 冰川 流速 的 时 间 变 化 ,本 研究 沿 冰 
川中 心 线 提取 速度 剖面 图 (图 $)。 从 网 中 可 以 看 到 
冰川 中 上 游 在 2014 年 10 月 已 具有 较 高 的 流速 ,最 
高 流速 达 3.6 md"', 此 时 布 拉 尔 杜 冰川 处 于 高 速 运 
动 期 。 随 后 冰川 的 跃 动 前 锋 保持 较 高 的 流速 向 下 
游 推 进 , 在 此 过 程 中 最 高 流速 达 4.9 md'。 到 2015 
年 8 月 底 , 冰 川 表面 流速 急剧 下 降 , 跃 动 前 锋 的 推进 
随 着 冰川 流速 的 下 降 而 终止 ,冰川 进入 平静 期 。 然 
而 到 2016 年 1 月 初 ,冰川 表面 流速 突然 增加 ,最 高 
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Fig. 4 Elevation change 
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Fig. 5 Variation of daily average velocity along the mainstream of the Braldu Glacier 


流速 达到 3.9 md"', 这 种 高 流速 的 状态 维持 了 一 个 
H ,随后 冰川 进入 平静 期 。 结 合 ITS_LIVE 项 目 
(https://nsidc.org/apps/itslive/) 中 的 高 亚洲 年 均 流 速 
数据 对 布 拉 尔 杜 冰川 进行 分 析 ( 图 6) ,发 现 从 2006 
年 起 , 布 拉 尔 杜 冰川 主干 中 上 游 的 流速 逐年 升 高 ， 
有 具 有 明显 的 跃 动 前 锋 并 向 下 游 缓慢 推进 ,此 时 冰川 
已 经 显示 出 跃 动 特征 趋势 。 布 拉 尔 杜 冰 川 在 
2012—2015 年 处 于 快速 运动 期 ,在 此 期 间 最 高 流速 
可 达 400 ma'。 布 拉 尔 杜 冰川 由 多 个 支 冰 川 汇流 
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图 6 1989 一 2018 年 布 拉 尔 杜 冰 川 沿 主流 线 流速 变化 
Fig.6 Velocity variation along the mainstream of the Braldu 
Glacier during 1989—2018 


而 成 ,图 7 为 布 拉 和 尔 杜 冰川 详细 的 冰川 流速 数据 ,从 
图 7a~d、 图 7g 中 可 以 观察 到 在 跃 动 期 间 不 仅 冰 川 主 
干 具有 较 高 的 流速 , 支 冰川 流速 也 相对 较 高 ;而 在 
平静 期 , 布 拉 尔 杜 冰川 主干 和 其 支 冰川 均 处 于 缓慢 
流速 的 状态 (图 7e~f、 图 7 h~i)。 
3.2 冰川 表面 高 程 变化 

为 了 定量 分 析 冰 川 跃 动 前 后 表面 高 程 的 变化 ， 
本 文 基于 大 地 测量 法 获得 2000 一 2014 年 .2014 一 
2018 年 布 拉 尔 杜 冰川 的 表面 高 程 变 化 (图 8), 由 于 
2014 年 的 TSX/TDX 数据 没有 完全 覆盖 布 拉 尔 杜 冰 
川 ,因此 导致 图 8 部 分 区 域 缺 失 。2000 一 2014 年 3、 
4 号 支 冰川 向 主干 输送 大 量 的 冰川 物质 ,人 致 使 冰川 
主干 中 上 游 有 明显 隆起 ,1、2、5、6 号 支 冰川 与 主干 
相 接 的 位 置 也 明显 增 厚 , 而 冰川 主干 的 下 游 以 消融 
为 主 ,其 高 程 明显 降低 。2014 一 2018 年 布 拉 尔 杜 冰 
川 主 干 中 上 游 以 及 各 支 冰 川 有 不 同 程 度 的 减 薄 , 而 
冰川 主干 的 接收 区 明显 增 厚 ,最 大 增 厚 可 达 120 mo 

沿 冰 川中 心 线 提取 布 拉 尔 杜 冰川 2000 一 2014 
年 和 2014 一 2018 年 的 表面 高 程 变 化 (图 9) ,结果 表 
明 在 2000 一 2014 年 ,冰川 的 下 游 (0~13 km) 以 消融 
为 主 ,最 大 减 薄 89 m, 而 冰川 主干 中 上 游 ( 距 末端 
13~26 km) 有 不 同 程度 增 厚 ;2014 一 2018 年 ,冰川 主 
干 的 接收 区 (6~13 kmy) 高 程 显著 增加 ,冰川 主干 中 
上 游 减 薄 5~20 m 不 等 。 
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图 7 2014 一 2017 年 布 拉 尔 杜 冰川 详细 的 流速 变化 
Fig.7 Detailed velocity changes of the Braldu Glacier during 2014—2017 


I 6 号 支 冰川 从 2006 年 起 便 处 于 跃 动 状态 ,这 与 
4 讨论 Gardelle 等 中 研究 的 结果 相 一 致 ,但 他 们 的 研究 中 
ITS_LIVE 年 均 流速 (图 6) 表 明 此 冰川 从 2006 认为 1 号 支 冰 川 也 处 在 跃 动 期 ,而 本 文 仅 通 过 


年 开始 缓慢 加 速 ,并 逐渐 向 下 游 推进 。 根 据 Landsat 2000 一 2014 年 高 程 和 表 碘 的 变化 无 法 判别 1 号 支 冰 
影像 发 现 5.6 号 支 冰川 的 冰 碎 从 2006 年 起 逐渐 向 。 川 是 否 从 2006 年 开始 跃 动 。 从 2009 年 开始 ,3、4 号 
外 扩展 ,结合 图 8a 中 的 高 程 变 化 可 以 发 现 5.6 号 支 SCM BAR ig 23 km 处) 逐渐 扩大 并 沿 主干 
冰川 与 冰川 主干 相 接 位 置 有 明显 的 增 厚 , 故 认为 5、 向 下 游 推 进 , 结 合 冰 川 表面 流速 的 结果 (图 5、 图 6) 
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图 8 布 拉 尔 杜 冰川 表面 高 程 变化 
Fig. 8 Change of the Braldu Glacier surface elevation 
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图 9 布 拉 尔 杜 冰川 沿 主流 线 表 面 高 程 及 高 程 变 化 


Fig.9 Elevation and elevation change profile of the Braldu 


Glacier along the mainstream 


发 现 布 拉 尔 杜 冰川 从 2011 年 开始 处 于 高 速 运动 的 
状态 , 随 着 支 冰 川 和 主干 上 游 的 冰川 物质 向 下 游 推 
进 , 冰 川 主干 中 部 的 表 磺 受到 挤 压 发 生 扭 曲 并 向 下 
游 继续 推进 。 而 布 拉 尔 村 冰川 表面 冰 厂 的 变化 是 
由 于 表面 速度 的 不 均匀 性 所 导致 ,同时 也 表明 冰川 
SbF ERSTE’ 

HEE eh AK J] ee TE Bat E A AKE KS ZEAE eB A FR 


杜 各 支 冰 川 跃 动 的 起 始 时 间 并 不 一 致 ,根据 2000 一 
2014 年 布 拉 尔 杜 冰川 高 程 变化 结果 (图 8a) 可 以 推 
断 布 拉 和 尔 杜 冰川 的 各 分 支 在 2014 年 之 前 均 已 开始 
跃 动 。 图 8a 的 结果 显示 在 2014 年 之 前 仅 有 3、4 号 
支 冰川 的 物质 向 下 游 传 递 ,化 使 冰川 主干 的 中 上 游 
(13~20 km) 明 显 增 厚 ;1、2、5、6 号 支 冰川 虽 在 此 之 
前 已 经 发 生 跃 动 ,但 来 自 这 四 条 支 冰川 上 游 的 物质 
大 多 聚集 在 与 冰川 主干 相 接 的 位 置 ,故此 认为 3、4 
号 支 冰川 的 跃 动 引 起 主干 冰川 的 第 一 阶段 跃 动 。 
而 在 2014 年 之 后 ,冰川 主干 的 中 上 游 和 各 支 冰川 的 
冰川 物质 沿 冰川 主干 向 下 游 传播 ,导致 冰川 主干 的 
接收 区 (6~13 km ) 的 高 程 显著 增加 ,最 大 增 厚 约 120 
m, 而 1.2.5.6 号 支 冰川 明显 减 薄 , 由 此 认为 1.2.5、 
6 号 支 冰川 的 跃 动 引起 第 二 阶段 的 跃 动 。 

由 于 缺乏 2006 一 2013 年 详细 的 流速 数据 ,无 法 
HR UB Fi F ZR FE DK JII BR 5h W H tak E 2E AE BF IE o 
Quincey 等 ”利用 光学 影像 特征 匹配 法 获取 2014 年 
布 拉 尔 杜 冰川 夏季 表面 流速 数据 ,发 现 8 月 的 流速 
与 7 月 相 比 大 幅 下 降 , 且 通过 现 有 的 流速 数据 也 发 
现 此 冰川 在 2015 年 8 月 底 突然 减速 ,在 保持 一 段 低 
速 运动 后 ,于 2015 一 2016 年 冬季 突然 加 速 , 随 后 逐 
渐 减 速 。 布 拉 尔 村 冰川 流速 的 季节 性 变化 比较 符 
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合 冰 下 水 文 控 制 冰 川 的 跃 动 特征 ,速度 在 夏季 急剧 
下 降 , 在 冬季 会 上 升 下 。 受 水 文 控 制 的 冰川 跃 动 ， 
流速 与 水 下 水 压 直 接 相关 ,流速 的 大 小 取决 于 冰 下 
排水 系统 的 状态 5 。 在 夏季 ,大 量 融 水 的 输入 致使 
排水 渠道 高 效 的 运转 , 冰 下 水 压 下 降 ,导致 流速 下 
降 。 由 于 冰 体 的 持续 压缩 运动 和 可 塑性 ,排水 通道 
在 秋冬 季节 随 着 融 水 的 减少 而 变 得 低 效 运转 ,随后 
冰 下 水 压 升 高 ,最 终 冬 季 冰 流 加 速 。 

与 喀 喇 昆仑 其 他 跃 动 冰川 相 比 , 布 拉 尔 杜 冰川 
在 跃 动 期 间 并 未 表现 出 较 高 的 流速 ,最 高 流速 仅 
4.9 md"'。 而 与 布 拉 尔 杜 冰川 相 邻 的 Hispar 冰川 在 
2015 年 春季 的 最 高 流速 可 达 14 md', 此 冰川 在 
2015 年 夏季 被 观察 到 流速 突然 下 降 , 却 在 2015 一 
2016 年 冬季 再 次 加 速 ,其 跃 动 也 被 认为 是 与 冰 下 水 
文 状况 有 关 , H. Hispar 冰 川 此 次 跃 动 主要 受 其 支流 
B, Quincey 等 ”通过 对 一 些 喀 喇 昆仑 冰川 跃 
动 过 程 研究 发 现 , 喀 喇 昆仑 地 区 的 冰川 跃 动 并 不 完 
全 符合 经 典 的 水 文 控制 模型 或 热 控制 模型 ,地形 和 
表 奢 物 也 是 影响 跃 动 的 重要 因素 。 布 拉 尔 杜 冰川 
表面 覆盖 着 厚 厚 的 表 磺 和 和 死 冰 ,这 些 物质 可 能 会 阻 
碍 冰川 物质 向 下 游 输 送 。 

此 外 , 张 祥 松 曾 观察 到 布 拉 尔 村 冰川 的 末端 
在 1968 一 1973 年 前 进 了 约 900 m, JF ERTZE YKT 
与 书 斯 米 亚 兹 河口 的 冲积 局 相 接 ;Rankl 等 ”利用 
SRTM TerraSAR-X .ALOS PALSAR 和 ENVISAT ASAR 
等 多 源 雷 达 影 像 对 喀 喇 昆仑 冰川 变化 进行 研究 ,发 
现 布 拉 尔 杜 冰川 上 次 跃 动 发 生 在 1976 年 之 前 ,因此 
推断 布 拉 尔 杜 冰川 的 跃 动 间 隔 约 为 40 ao 


5 结论 

(1) 根据 ITS_LIVE 年 均 流 速 数据 发 现 ,从 2006 
年 开始 , 布 拉 尔 杜 冰川 年 均 流 速 逐渐 增加 ,展现 出 
跃 动 趋势 。 

(2) 结合 布 拉 尔 杜 冰川 高 程 变化 数据 发 现 , 在 
2014 年 之 前 ,此 冰川 的 各 支 冰 川 均 已 开始 路 动 , 但 
仅 有 3、4 号 支 冰川 持续 向 下 游 进 行 物 质 传 输 , 而 其 
余 支 冰川 在 2014 年 之 后 开始 向 下 游 进 行 物质 的 传 
输 , 可 认为 由 于 3、4 号 支 冰川 的 跃 动 引 起 主干 冰川 
的 第 一 阶段 跃 动 ,1、2、5、6 分 支 在 2014 一 2018 年 物 
质 向 下 游 迁 移 引 起 第 二 阶段 跃 动 。 

(3) 根据 布 拉 尔 杜 冰川 2014 一 2018 年 详细 的 
流速 数据 ,发 现 此 冰川 流速 的 变化 特征 比较 符合 水 
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文 控 制 冰 川 的 跃 动 特征 。 结 合 现 有 的 文献 及 影像 
数据 ,推断 布 拉 尔 杜 冰 川 的 跃 动 间 隔 约 为 40 ao 
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Characteristics of recent surging of Braldu Glacier, Central Karakoram 


WANG Zhenfeng'’, JIANG Zongli’, LIU Shiyin’, 

YANG Jingrui, MA Zhiyuan’, WANG Wenging’ 
(1. School of Resource & Environment and Safety Engineering, Hunan University of Science and Technology, Xiangtan 411201, 
Hunan, China; 2. School of Earth Sciences and Spatial Information Engineering, Hunan University of Science and Technology, 


Xiangtan 411201, Hunan, China; 3. Institute of International Rivers and Eco-Security, Yunnan University, 
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Abstract: The Karakoram hosts a large number of surge-type glaciers. However, the mechanisms of glacier surges 
are not clear yet. In this paper, Landsat, Sentinel-1A, TSX/TDX, and other multisource remote sensing data were 
used to obtain the changes in surface elevation and surface velocity of the Braldu Glacier during the surge period. 
The results show that the velocity of the Braldu Glacier has gradually increased since 2006. After 2013, it entered 
a period of rapid movement, with the highest velocity reaching 4.9 m:d '. At the end of August 2015, the surface 
velocity decreased sharply. The glacier then kept a low velocity until January 2016, when the velocity increased 
again. The rapid movement lasted only for about a month and then entered a quiet period. From 2000 to 2014, 
there was a significant uplift in the middle and upper reaches of the glacier trunk, whereas the glacier receiving 
area decreased significantly, with a maximum thinning of 89m. From 2014 to 2018, the upper and middle 
reaches of the glacier trunk and its tributaries were thinned to varying degrees, and the elevation of the receiving 
area of the glacier trunk increased significantly, with a maximum thickness increase of 120 m. According to the 
characteristics of surface elevation change and velocity change, it is considered that the tributary of the Braldu 
Glacier caused this surge, which in turn is greatly affected by hydrological mechanisms. Combined with available 
data and literature, it is inferred that the interval of the Braldu Glacier’s surge is about 40 years. 


Key words: glacier surging; surface velocity; surface elevation; TanDEM-X; Braldu Glacier 


